First- and Second-Order Phase Transitions, Fulde-Ferrel Inhomogeneous State and 
Quantum Criticality in Ferromagnet/Superconductor Double Tunnel Junctions 
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First- and second-order phase transitions, Fulde-Ferrel (FF) inhomogeneous superconducting (SC) 
state and quantum criticality in ferromagnet/superconductor/ferromagnet double tunnel junctions 
are investigated. For the antiparallel alignment of magnetizations, it is shown that a first-order phase 
transition from the homogeneous BCS state to the inhomogeneous FF state occurs at a certain bias 
voltage V* ; while the transitions from the BCS state and the FF state to the normal state at V c are 
of the second-order. A phase diagram for the central superconductor is presented. In addition, a 
quantum critical point (QCP), Vqcp, is identified. It is uncovered that near the QCP, the SC gap, 
the chemical potential shift induced by the spin accumulation, and the difference of free energies 
between the SC and normal states vanish as |V — Vqcp\ zv with the quantum critical exponents 
zv — 1/2, 1 and 2, respectively. The tunnel conductance and magnetoresistance are also discussed. 

PACS numbers: 73.40.Gk, 75.70.Pa, 73.40.Rw 



Introduction. — Spin-dependent transport plays an 
essential role in magnetic hybrid nanostructures in 
the field of spintronics (see, e.g. Refs. 0, 0, 
II 3 S II II for review). Among others, ferro- 
magnet/superconductor (F/S) heterostructures have at- 
tracted much attention theoretically |il fToL [III fl2l H^. 
Q and experimentallydl H [H M, (H liflW^ 1 re- 
cent years. For F/S/F double tunnel junctions, it has 
been observed that the superconductivity is suppressed 
by the injection of spin-polarized current (e.g. Refs. 
SmElUQElIi), that is due to the nonequilib- 
rium spin accumulation. When spin-polarized electrons 
are injected into the superconductor, a spin density is 
accumulated near the interfaces owing to the spin im- 
balance, thereby giving rise to an equivalent, small mag- 
netic field that acts as a pair-breaking field, which leads 
to a suppression of superconductivity. There has been a 
recent studyQ showing that the homogeneous supercon- 
ducting (SC) state is strongly suppressed with increasing 
the bias voltage and completely destroyed at a critical 
voltage by the nonequilibrium spin density in the antipar- 
allel alignment of magnetizations. This study is qualita- 
tively consistent with the experimental observation for 
high biases, but is inconsistent for low biases 18J. 

On the other hand, about forty years ago, Fulde and 
Ferrel (FF)[3], and Larkin and Ovchinnikov (LO)[ij| 
independently, found that the SC order parameter can 
be modulated in real space by a spin-exchange field of 
a ferromagnet. Later, such an inhomogeneous supercon- 
ducting state has been extensively explored under various 
circumstances (e.g. [l^.l24l|'). As the nonequilibrium spin 
accumulation may lead to an equivalent magnetic field 
in the central supercondutor, the FFLO state, which is 
simply omitted in the previous treatment 9], might be 
inevitable in the F/S/F double tunnel junction. 

To understand profoundly the spin-dependent trans- 
port properties, in this Letter, the F/S/F double tunnel 
junction shall be systematically revisited. It is shown 



that in the antiparallel aligment, a first-order phase tran- 
sition from a homogeneous Bardeen-Cooper-Schricffcr 
(BCS) SC state to the inhomogeneous FF SC state oc- 
curs at a certain bias voltage, while the transitions from 
the BCS state and the FF state to the normal state are of 
the second-order. A phase diagram for the central super- 
conductor is identified. Besides, a quantum critical point 
(QCP) is specified at bias Vqcp, near which the SC or- 
der parameter, the chemical potential shift induced by 
the spin accumulation, and the difference of free energies 
between the SC and normal states vanish as \V—Vqcp\ zl ' 
with the quantum critical exponents zv = 1/2, 1 and 2, 
respectively. The tunnel conductance and magnetoresis- 
tance are also obtained. 

Model. — Consider a symmetric F/S/F double tun- 
nel junction with the left and right ferromagnetic (FM) 
electrodes applied by bias voltages —V/2 and V/2, re- 
spectively. The two identical FM electrodes are sepa- 
rated from the central superconductor by two insulating 
thin films. The central superconductor is presumably 
described within the framework of BCS theory. Sup- 
pose that the energy relaxation time of quasiparticles 
is shorter than the tunneling time, while the latter is 
shorter than the spin relaxation time. As the resistance 
of this tunnel junction with insulating thin films is greater 
than that of a conventional metallic contact, the Andreev 
reflection effect can be reasonably ignored for simplicity. 

From the standard tunneling Hamiltonian, and in 
light of the linear response theory, the tunneling current 
through the jth junction can be readily obtained by 
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where T is the tunneling matrix element, j = 1, 2, Dj a is 
the subband density of states (DOS) in the jth FM elec- 
trode, cr = ±1 for spin up and down, respectively, and 
rji = 1, 772 = — 1. The quantities S, Q, N and N are 
defined by 
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where /o(z) denotes the Fermi distribution function of 
thermal equilibrium in FM electrodes, /k CT is the nonequi- 
librium distribution function of quasiparticles with en- 
ergy i?k and spin <j (=Ti I) m the central superconductor, 
S and Q represent the spin density and the quasiparticle 
charge density, describing the spin imbalance and quasi- 
particle charge imbalance in the central superconductor, 
respectively. 

Fulde-Ferrel State. — As there appears the nonequilib- 
rium spin accumulation near the interfaces of the tunnel 
junction, it is believed that the SC state would gener- 
ally include the inhomogeneous FFLO state |22j, |23| in 
addition to the homogeneous BCS state. Without loss of 
generality, we suppose that the SC order parameter takes 
the form of FF type 22]: A(r) = A q e iqr with A q the am- 
plitude of the order parameter. When q = 0, it recovers 
the homogeneous BCS state. The quasiparticle disper- 
sion i?k, the coherence factors iik and are given by 



and vl = |(1 - £ k / J(£ + A 2 ), respectively, where £ k is 
the free electron energy relative to the chemical potential, 
vf the Fermi velocity and x = k • q/ (qk) the cosine of 
the angle between k and momentum q of a Cooper pair. 
The amplitude A q is determined by the gap equation 
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where Vbcs is the BCS-type pair interaction. The values 
of q will be specified later. 

Let us proceed to determine the nonequilibrium distri- 
bution function ■ In the absence of spin-flip scatter- 
ing, the spin up and down tunneling currents are indepen- 
dent, and should be conserved, i.e. l\ a — l2 a , yielding 



S F = for the parallel alignment, (7) 
S A = PN A for the antiparallel alignment, (8) 
Q F = Q A = for both alignments, (9) 

where the superscripts F and A refer to the paral- 
lel and antiparallel alignments, respectively, and P = 
— Dji \ /(Dji + Dji) is the spin polarization of the 



FM electrodes. These solutions show that the nonequi- 
librium spin accumulation exists only in the antipar- 
allel configuration. In the above derivation, we have 
adopted the conventional constant DOS approximation, 
£ k (- ' ' ) - N(0) J™° d (• • • )d£jk, where w D is the cut-off 
(Debye) energy, and iV(0) denotes the DOS of free elec- 
trons at the Fermi level. The quantity N vanishes iden- 
tically for both alignments of magnetizations since the 
integrand is an odd function of £k- Eq. © requires that 
should differ from in the presence of the tunnel- 
ing current. Following Ref.jtJ, we consider the solutions 
of the form: 
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where Sfj, is introduced as the chemical potential shift 
induced by the nonequilibrium spin accumulation, and 
plays essentially the same role as the spin-exchange field 
explored by FF in their seminal article!^. This kind of 
solutions may be applicable if the thickness of the central 
superconductor is much smaller than the spin diffusion 
length, and the spin relaxation time is sufficiently long. 
Eqs. ((THf and (fTTJ) are the solutions of Eqs. and 
©. However, Eqs. $ and © with Eqs. should be 
solved in a self-consistent manner to specify <5/i and A q as 
functions of the bias V, temperature T and polarization 
P. 

For these coupled equations, when the self-consistent 
multiple solutions corresponding to different q appear, 
only the value of q that leads to the lowest free energy of 
the system is retained. In an inhomogeneous supercon- 
ducting state, the nonzero solution for A q implies only 
the local minimum of the free energy, which does not nec- 
essarily mean the stable state. In order to clarify this 
issue, one must compare the free energies of the homo- 
geneous BCS, the inhomogeneous FF and normal states, 
as emphasized by Abrikosov|2^. The free energy of the 
present system can be obtained by integrating the gap 
eciuation[2fij 
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where Fg^&nd F^ A ^ stand for the free energy of the 
SC state and the N state, respectively, and S/j, f = 0, 
5fi A = S[i. It turns out that both the homogeneous BCS 
(q = 0) and the FF (q ^ 0) SC solutions are possible 
in the antiparallel alignment, while the former solution 
is always favorable in the parallel configuration. 
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Results. — In the parallel alignment of magnetizations, 
since there is no spin and charge accumulation in this 
circumstance, the SC order parameter does not depend 
on the bias voltage; while in the antiparallel configura- 
tion, the situation becomes complicated, as the nonequi- 
librium spin accumulation characterized by S A intervenes 
in. Figure 1 presents the bias voltage dependence of the 
SC order parameter and the chemical potential shift at 
T/T c = 0.2, with T c the SC critical temperature, for the 
antiparallel alignment. It is observed that the order pa- 
rameter A A remains almost constant at low biases and 
is in the homogeneous BCS state till a specific bias volt- 
age V* = 1.36Ao/e at which A A drops suddenly, where 
P = 0.4, and Ao is the BCS zero-temperature energy 
gap. Then, A A goes into the inhomogeneous FF state 
with q ^ 0, decreases with the bias, and vanishes com- 
pletely at V — V c where superconductivity is quenched, 

1. e. A q (V c ) = 0, as shown in Fig. 1(a). At V = V*, 
there is a discontinuity for A A , implying that a first- 
order phase transition from the homogeneous BCS phase 
to the FF phase exists in the system. The present result 
is quite different from that given in Ref.Q where the in- 
homogeneous FF state was simply ignored. The chemical 
potential shift grows in the homogeneous BCS state with 
increasing the bias, and exbihits a jump at V* , then in- 
creases slowly in the FF state till V c , and is linear for 

V ^ V c in the normal state, as shown in Fig. 1(b). We 
have found that the free energy of the inhomogeneous 
FF phase is lower than the homogeneous BCS phase for 
V c > V ^ V*, suggesting that the FF state is stable. At 

V = V* , the free energy of the FF state coincides with 
that of the homogeneous BCS state, revealing that the 
FF state can coexist with the BCS state at V*. 

The magnitude of momentum q of a Cooper pair is 
in general a function of bias voltage, as shown in Fig. 

2. At V < V* , the system is in the homogeneous BCS 
state, and thus |q| is zero; while for V* ^ V < V c , |q| 
corresponding to the lowest free energy varies nonmono- 
tonically with the bias, implying that the momenta of 
Cooper pairs in the stable, inhomogeneous SC state are 
not fixed. Note that there is a discontinuity for |q| at 

V = V* 1 which is again the signature of the first-order 
phase transition. 

Phase Diagram. — In the antiparallel configuration, a 
schematic phase diagram for the central superconductor 
in the T — V plane could be obtained, as shown in Fig. 3 
for P = 0.4. There exist three phases: the homogeneous 
BCS phase with A constant (q = 0) at low bias; the 
inhomogeneous FF phase with q ^ for V* ^ V < V c 
at low temperature; the normal phase for V ^ V c . Along 
the phase boundary from B to C, the homogeneous BCS 
phase coexists with the spatially modulated FF phase; 
along the V c boundary where the SC gap closes, there are 
two transitions: one is from the homogeneous BCS state 
to the normal state, and the other is from the spatially 
modulated FF state to the normal state, which are of 



1.0 



< 

^ 0.5 



0.0 



0.6 



<° 0.3 



0.0 





to 

\. FF 


BCS I 


P=0.4 




T/T =0.2 

' C 






(b) 


FF 


- BCS/ 


' N P=0.4 

T/T =0.2 



0.0 0.5 1.0 1.5 

eV/2A g 

FIG. 1: (Color online) The bias dependence of the order pa- 
rameter A A (a) and the chemical potential shift (b) in the an- 
tiparallel alignment of magnetizations at P = 0.4, T/T c = 0.2, 
where a discontinuity from the first-order phase transition is 
observed. The dashed line is for the normal state. 




FIG. 2: (Color online) The bias dependence of the momen- 
tum of a Cooper pair in the antiparallel alignment of mag- 
netizations at P = 0.4, T/T c = 0.2, where a discontinuity is 
observed. 



second-order. At point C, three phases meet, implying it 
can be viewed as a Lifshitz voint |27|: at point QCP where 
V QC p satisfies PeV QC p/2A = 0.754 [= S(j.(Vqcp)], a 
quantum phase transition (QPT) from the ordered FF 
state to the disordered normal state occurs at T = 0. 
TMR. — The total currents are given by I F = I$N F ', 

2 



T 



and I A = I Q (1-P 2 )N A , where I = 2ne 

and A^" 4 ) is given by Eq. Q. The differential conduc- 
tance G FlyA ^ and the tunneling magneoresistance TMR 
are obtained by G F ( A ' = 



dV 



and TMR = g, 



1. 



The bias dependence of the tunnel conductance and the 
TMR is shown in Fig. 4. The behavior of G F is consis- 
tent with that of Ref.Q, but G A and TMR differ from 
those in Ref. owing to the intervention of the FF state, 
where the oscillating behaviors for G A and TMR are 
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FIG. 3: (Color online) A schematic phase diagram of the 
central superconductor in the T-V plane for the antiparallel 
alignment of magnetizations at P = 0.4. Three phases are 
observed: the homogeneous BCS state is separated by the 
V* boundary line via a first-order phase transition from the 
inhomogeneous FF state, while the SC states are separated by 
the V c boundary line via a second-order phase transition from 
the normal state. A quantum phase transition is sepcified at 
Vqcp (see context). 



As V approaches Vqcp, the amplitude of SC order pa- 
rameter, the chemical potential shift induced by the spin 
accumulation, as well as the difference of free energies 
between the SC and normal states vanish as A q /A , 
[SKVqcp) ~ Sfi(V)], f£~F$~\V- V QCP r with the 
critical exponents zv = 1/2, 1 and 2. Obviously, such a 
QPT is of second-order. The present system offers a nice 
example of the QPT. 

In summary, we have revisited the spin-dependent 
transport in F/S/F double tunnel junctions by taking the 
proximity effect into account. It is found that in the an- 
tiparallel configuration, the first- and second-order phase 
transitions, the inhomogeneous FF state, and quantum 
criticality can be revealed simultaneously in the central 
superconductor. The present study rectifies the previous 
result^ where the FF state was simply ignored. 

This work is supported in part by the NSFC (Grant 
Nos. 90103023, 10104015, 10247002). 
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FIG. 4: (Color online) The bias dependence of the tunnel con- 
danctance G F ' A (a) and the magnetoresistance (b) at P = 0.4, 
T/T c = 0.2, where is the conductance in the normal state 
for the antiparallel alignment. The oscillation characteristics 
are observed. 



seen. These oscillation characteristics could be served as 
good tests for experimentally observing the inhomoge- 
neous FF state, and are qualitatively consistent with the 
recent observation 20J . 

Quantum Criticality. — Since a phase transition from 
the inhomogeneous FF state to the normal state appears 
at Vqcp, which just happens at T = 0, it is nothing 
but a manifestation of the QPTj2a|. At V < Vqcp, the 
system is in an ordered ground state; while V > Vqcp, 
it is in a disordered ground state, suggesting Vqcp is 
the quantum critical point, where the SC gap vanishes. 
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